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Observations of nearby galaxies have firmly established, over a broad range of galactic environments 
and metallicities, that star formation occurs exclusively in the molecular phase of the interstellar 
medium (ISM). Theoretical models show that this association results from the correlation between 
chemical phase, shielding, and temperature. Interstellar gas converts from atomic to molecular only 
in regions that are well shielded from interstellar ultraviolet (UV) photons, and since UV photons are 
also the dominant source of interstellar heating, only in these shielded regions does the gas become cold 
enough to be subject to Jeans instability. However, while the equilibrium temperature and chemical 
state of interstellar gas are well-correlated, the time scale required to reach chemical equilibrium is 
much longer than that required to reach thermal equilibrium, and both timescales are metallicity- 
dependent. Here I show that the difference in time scales implies that, at metallicities below a 
few percent of the Solar value, well-shielded gas will reach low temperatures and proceed to star 
formation before the bulk of it is able to convert from atomic to molecular. As a result, at extremely 
low metallicities, star formation will occur in a cold atomic phase of the ISM rather than a molecular 
phase. I calculate the observable consequences of this result for star formation in low metallicity 
galaxies, and I discuss how some current numerical models for H2-regulated star-formation may need 
to be modified. 

Subject headings: galaxies: star formation — ISM: atoms — ISM: clouds — ISM: molecules — stars: 
formation 



1. INTRODUCTION 

In present day galaxies, star formation is very well- 
correlated with th e molecular phase of the interstel- 
lar medium QSM) rtWong fc Blitdl2002t iKennicutt et al.l 
120071 : iLerov et al.ll2008t iBigiel et all 120081 ). In contrast, 
in the inner parts of disks where there are significant 
molecular fractions, star formation correlates very poorly 
or not at all with the atomic ISM. At large galctocen- 
tric radii where the ISM becomes atomic-dominated star 
formation does begin to correlate with H I, but this ap- 
pears to be only because H2 itself becomes correlated 
with H 1, and the H2 forms stars in the same way regard- 
less of where it is found within a galaxy (Bigi el et al.l 
l2010t iSchruba et al.l 120111) . Strong association between 
star formation and H2 and a lack of association with H I is 
also found down the lowest metallicity syst ems that have 
been measured, at roughly 20% of Solar (jBolatto et al.l 
2011). In summary, all available observational data indi- 
cates that star formation occurs only where the hydrogen 
in the ISM has converted to H2. 

Theoretical models have explained these observa- 
tions as resulting f rom a correla t ion between chemistry 
and temperature (ISchavd 120041: IKrumholz et al.l 120111: 
IGlover fe Clark! l2012aD ~Molecular hydro gen is not an 
important coolant in modern-day galaxies, and while car- 
bon monoxide (which forms only when it is catalyzed by 
H 2 - Ivan Dishoeck fc BlackHT9 86) is, the C 11 found in H I 
regions is almost as effective. However, H2 is an excellent 
proxy for the presence of cold gas because both are sen- 
sitive to destruction by UV photons, which photodisso- 
ciate H2 and increase the temperature through the grain 
photoelectric effect. As a result, both H2 and low tem- 
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perature gas are found only in regions of high extinction 
where the UV photon density is far below its mean value 
in the ISM, and, conversely, any region that where the 
photodissociation rate is high enough to convert the bulk 
of the ISM to H 1 is also likely to be warm. Since low 
temperatures that remove thermal pressure support are 
a prerequisite for collapse into stars, this correlation be- 
tween temperature and chemical state in turn induces a 
correlation between star formation and chemical state. 

However, the correlation between H2 and star forma- 
tion must break down at sufficiently low metallicities. 
Before the first stars formed in the universe, and for a 
short time thereafter, there were no or very few heavy 
elements. As a result, forming H2 was extremely dif- 
ficult due to a lack of dust grain surfaces to catalyze 
the H 1 — > H2 reaction. Theoretical models of star for- 
mation in such environments indicate that H2 fractions 
remain extremely small until the density rises so high 

^ 10 9 cm~ 3 ) tha t H2 can form v i a three-body reactions 
Palla et al.1 Il98l ILepp . fc Shulll fl98i lAhn fc Shapirol 
l2007tlOmukai et al.ll2010D . Theunderlving physical basis 
for this result is a disconnect of timescales: the equilib- 
rium chemical state the gas would reach after a very long 
time would be ^-dominated, but the cooling and star 
formation times are short enough that the gas does not 
reach equilibrium before collapsing into a star. 

While this result has been known for zero and ex- 
tremely low metallicity systems for some time, the re- 
lationship between chemical state and star formation in 
intermediate metallicity regime, for which observations 
are at least in principle possibl e in the local univers e, has 
received fairly little attention. lOmukai et al.l (|2010T ) con- 
sider the chemical evolution of collapsing gas cores with 
metallicities from to Solar, and investigate under what 
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circumstances they can form H2. However, because their 
calculation starts with unstable, collapsing cores, it does 
not address the question of in what phase of the ISM one 
expects to find such collapsing regions in the first place, 
which is the central problem for understanding the ob- 
served galactic-scale corr elation between ISM ch emical 
state and star formation. iGlover fc Clarkl (|2012bD simu- 
late the non-equilibrium chemical and thermal behavior 
of clouds with metallicities from 1% of Solar to Solar. 
They find that the bulk of the cloud material converts to 
H2 before star formation in the high metallicity clouds 
but not in the lowest metallicity ones, indicating that 
the star formation - H2 correlation should begin to break 
down at metallicities observable in nearby galaxies. How- 
ever, given the computational cost of their simulations, 
they are able to explore a very limited number of cases, 
and it is unclear how general their results might be. 

The goal of this pap er is to go beyond t he stu dies 
of lOmukai et~aT1 (|2010h and IGlover fc Clarkl (|2012bl) by 
deriving general results about the correlation between 
chemical state and star formation over a wide range of 
environments and metallicities. I do not perform de- 
tailed simulatio ns, such as those of lOmukai et al.l and 
IGlover fc ClaTkL for every case. Instead, I rely on fairly 
simple models that can be integrated semi-analytically. 
The benefit of this approach is that it is the only way to 
survey a large parameter space, and thereby to answer 
the central questions with which I am concerned: under 
what conditions do we expect the correlation between 
star formation and H 2 to break down? When such a 
breakdown occurs, what is the governing physical mecha- 
nism that causes it? What are the resulting observational 
signatures? What are the implications of this breakdown 
for the models of star formation commonly adopted in 
studies of galaxy formation? In the remainder of this 
paper, I seek to answer these questions. 

2. MODEL 

Consider spatially uniform gas characterized by a mean 
number density of H nuclei t!h, column density of hy- 
drogen nuclei N-r, metallicity Z' relative to Solar, and 
temperature T. A fraction /h 2 of the H nuclei are 
locked in H2 molecules. It is generally more conve- 
nient to characterize models by values of the visual ex- 
tinction Ay instead of -/Vh- These two are related by 
Av/Nn ~ 4.0 x 10~ 22 Z' mag cm 2 , with the normal- 
ization chosen as a compromise between the values for 
Milky Way extinction and the extinction curves of the 
Large and Small Magellanic Clouds adjusted to Milky 
Way metallicityQ 

2.1. Timescale Estimates 

We are interested in following the behavior of initially 
warm, atomic gas, and considering whether it will be able 
to cool to temperatures low enough to allow star forma- 
tion, and how its chemical state will evolve as it does 
so. Before computing detailed evolutionary histories, it 
is helpful first to make a rough estimate of the timescales 

1 The value of Ay/Nn, a nd all other parameters in the following 
discussion, are taken from Krumholz ct al. (2011) unless stated 
otherwise; arguments for these choices are given in that paper. 
Ato mic data (Einstein c oefficients, collision ra tes) are all ta ken 
from Schoicr et al. (2005), and abundances from Drainc (2011). 



involved. In interstellar gas that is dense enough to be a 
candidate for star formation, but that is not yet molec- 
ular or forming stars, the dominant cooling process is 
emission in the [C 11] 158 /im line, which removes energy 
at a rate 

Acn ~ kcu-n5ckBTcuCn H (1) 

per H atom, where fccn-H ~ 8 x io~ 10 e -Tcn / T cm 3 s _1 
is the rate coefficient for collisional excitation of C 11 by 
H atoms, 5c ~ 1-1 X 10~ 4 Z' is the gas phase carbon 
abundance, Ten — 91 K is the energy of the excited C 11 
level over ks , and C = (n 2 ^) /n 2 ^ is a clumping factor that 
accounts for clumping of the medium on size scales be- 
low that on which we are computing the average. This 
expression assumes that C 11 collisional excitation is dom- 
inated by H rather than by free electrons, that the gas 
is optically thin, and that the density is far below the 
critical density for the line; I show below that all these 
assumptions are valid. The time required for the gas to 
reach thermal equilibrium is of order 
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Acn fc C n_H<5cXciiCnH 

= 0.036 f — ) e Tc "/ T 'Z'- 1 ^ 1 ^ 1 Myr, (2) 
\Tc11J 

where C\ = C/10 and no = nn/1 cm~ 3 . Note that the 
value of C will depend on the size scale over which the 
average density is defined; the fiducial value C = 10 is 
intermediate between the values C ~ 2 and C ~ 30 
that numerical experiments indicat e are best for ~ 10 
pc and ~ 100 pc scale, r espectively (jGnedin et alj|200l% 
IMac Low fc Gloveri[2012l) . 

Conversion of the gas from atomic to molecular form 
occurs primarily on the surface of dust g rains down to 
meta llicities as low as ~ 10~ 5 of Solar (Omukai et al. 
I2010D . This process occurs at a rate per H atom naTZC, 
where 1Z ss 3 x 10~ 17 Z' cm 3 s" 1 is the rate coefficie nt for 
H2 formation on grain surfaces (jWolfire et alJ l2008). The 
associated timescale for conversion of the gas to molecu- 
lar form is 
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= 105Z /_1 no l Cf 1 Myr. 



The ratio of the two timescales is 

^chem fcf!TT— HfcTnTT Ml 
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indicating that the gas will reach thermal equilibrium 
vastly before it reaches chemical equilibrium. 

This difference in timescale is only important if the 
cooling of gas toward thermal equilibrium is followed by 
star formation on a timescale that is too short for the 
conversion of atomic to molecular gas to keep up. It is 
therefore helpful to consider a third timescale: the free- 
fall time 



iff 



3tt 



32GnH/iHWH 



= 43n, 



-1/2 



Myr, 



(5) 



where /in ~ 1.4 is the mean mass per H nucleus in units 
of the hydrogen mass ran. This is the timescale over 
which star formation should begin once gas is gravita- 
tionally unstable. The time for which a cloud survives 
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after the onset of star formation is significantly uncer- 
tain, but even the longest modern esti mates are ~ lOtg , 
while some are as short as ~ ltn (lElmegreenl 120001: 
Tan et all 120061 : iKawamura et alJfeOOa iGoldbaum et al l 
201 lb . Comparing this timescale to the two previously 
computed gives 



mate form of lDraine fc Bertoldil ((19961 ) . 
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^chem 



8 .3 x 10-4 ( JL) e Tcn/T z ,-l c -l n -l/2 (g) 

\TciiJ 

(7) 



= 2.4Z / - 1 C 1 -V Va 



Thus we see that the thermal timescale will be smaller 
than the free-fall timescale down to extremely low metal- 
licities for reasonable ISM densities and temperatures, 
but the same cannot be said of the chemical timescale. 
For example, at n = 100, T = 1000 K, Z' = 10~ 3 , 
and C\ = 1, the above equation gives tthermAff = 1-0, 
while fchem/iff = 240. A cloud with these properties 
could cool and proceed to star formation on a free-fall 
timescale without difficulty, but would not build up a 
substantial amount of H2 until more than 100 free-fall 
times. If such a cloud were anything like the observed 
star-forming clouds in the Milky Way, it would likely 
have been destroyed by stellar feedback well before this 
point. Note that individual overdense regions within the 
cloud in the process of collapsing to stars would still con- 
vert to H2, since i c hemAff is a decreasing function of den- 
sity, dramatically so once three-body reactions begin to 
occur; however, since the star formation efficiency is low, 
the non-star-forming bulk of the cloud material would 
not, and thus the amount of H2 present per unit star 
formation would be greatly reduced. 

2.2. Chemical and Thermal Evolution Models 

The argument above is based on simple timescale esti- 
mates. To check whether the result is robust, it is neces- 
sary to construct more sophisticated cooling and chem- 
istry models. Below I describe a more detailed model, 
and how it may be evaluated numerically to follow the 
thermal and chemical behavior of a cloud. 

2.2.1. Chemical Evolution 

The H 1 to H2 transition is governed by two main pro- 
cesses: formation of H2 on the surfaces of dust grains 
and destruction of H2 by photodissociation. The former 
occurs at a rate per H atom nnlZ, where tt-h is the lo- 
cal (rather than average) number density and 1Z is the 
metallicity-dependent rate coefficient given in the main 
text. In a region of mean density nn , the mean rate per H 
nucleus at which H 1 converts to H2 is simply the number 
density-weighted average of the rate given above, which 
is Cn\CR-- The photodissociation rate per H2 molcule is 



diss 



Cdi SS ,oe- CT£iWH /shield C/Vh 2 ), 



(8) 



where Cdiss.o ~ 5 x 10 11 s 1 is the dissociation rate for 
unshielded gas (jDraine fc Bertoldil 1 1996T ). a d w 10~ 21 Z' 
cm~ 2 is the dust cross section per H nucleus for Lyman- 
Werner band photons, iVn is the column density of H 
nuclei, JVh 2 is the column density of H2 molecules, and 
/shield (^Vh 2 ) is the shielding function that describes H2 
self-shielding. For the latter quantity, I use the approxi- 
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0.035 



(l + 0.1a;/6 6 ) 2 (l + x) 0J 



,.5xl0~ 4 (l+:r) a 



( 9 ) 

where x — Nn 2 /5 x 10 14 cm 2 and be is the Doppler 
parameter for the gas in units of 10 6 cm s _1 ; I use 
&6 = 0.71, corresponding to a velocity dispersion of a 
5 km s _1 , roughly that observed in molecular clouds in 
nearby galaxies, but the results are quite insensitive to 
this choice. The value of Cdiss,o is that appropriate for the 
Milky Way's radiation field, and this value will vary from 
galaxy to galaxy and within galaxies. However, I show 
below that the dependence of the results on this choice is 
also quite weak, because of the exponential dependence 
of the dissociation rate on column density: a relatively 
large change in Co can be compe nsated for by a fa r smaller 
change in iVn or Nu 2 (also see lKrumholz et al"1l2011l ). 

Given these processes, the rate of change of the H 2 
fraction is given by 



dt 



/H 2 - (1 - fu 2 )nuCU - / H2 Cdiss(iVH, A^h 2 ). (10) 



Consistent with the simple uniform cloud assumption, 
I adopt 7Vh 2 = /h 2 ^Vh/2. Note that there is no factor 
of 2 or 1/2 in the second term due to a cancellation: 
there is one H 2 molecule per two H nuclei bound as H 2 
(multiplying by a factor of 1/2), but each dissociation 
generates two free H nuclei (multiplying by a factor of 
two). 

2.2.2. Thermal Evolution 

The thermal evolution depends on heating and cool- 
ing processes. For heating, the dominant mechanisms 
are cosmic ray heating and the grain photoelectric effect. 
The photoelectric heating rate per H nucleus is 



TpE — fpEn^'e 



(11) 



where Tpe.o ~ 4x 10~ 26 erg s _1 is the grain photoelectric 
heating rate in free space, and the numerical value is for a 
Milky Way radiation field. Note that, since photoelectric 
heating is dominated by photons with energies similar to 
the Lyman- Werner bands, the dust cross section ad here 
is the same as that used in the H 2 formation calculation. 
The heating rate per H nucleus from cosmic rays is 



Tcr — Ccr?cr 



(12) 



where Ccr is the primary cosmic ray ionization rate and 
qcr is the energy added per prima ry cosmic ray ioniza- 
tion For atomic gas, gcR ~ 6.5 eV (jDalgarno fc McCravl 
119721 ). The observed primary cosmic ray ionization 
rate in the Milky Way varies sharply between diffuse 
sightlines and dark clouds; the former show Ccr ~ 
3 x 10~ 16 s _1 , with roughly a dex dispersion, while 
those in dark clouds are an order of magnitude lower. 
Cm ^ 2 x 10~ 17 s- 1 dWolfire et alJ l2010t INeufeld et all 
120101: llndriolo fc McCalll I2012I ). Since cosmic ray heat- 
ing is only significant compared to photoelectric heat- 
ing in dark clouds, it seems more reasonable to adopt 
the latter as a fiducial value, although I verify below 
that this choice does not significantly affect the results. 
Observations also indicate that cosmic ray ionization 
rates vary roughly linearly with galactic star formation 
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rates (|Abdo et all [20Toh . For these reasons I adopt 
Ccr — 2 x 10~ 17 Z' s _1 as a fiducial value; the metal- 
licity scaling is a very rough way of accounting for the 
lower cosmic ray flux in galaxies with lower metalficities , 
masses, and star formation rates ()Krumholz et alj|2011l ). 

Cooling for interstellar gas with temperatures <C 10 4 
K is dominated by the 158 /im fine structure line of C II 
and the 63 and 145 /im fine structure lines of O I. The 
critical densities for these transitions are approximately 
4 x 10 3 , 6 x 10 4 , and 3 x 10 5 cm -3 , respectively assuming 
the dominant collision partner is H (see below). These 
critical densities are significantly higher than the densest 
cases I consider, and so I neglect collisional de-excitation 
compared to radiative de-excitation. Assuming all cloud 
atoms are in the ground state, the line-center optical 
depth of a cloud to photons emitted in one of these lines 
is 

r - —^ksxNn = [0.11, 0.38, 0.51}Z'A v . o b^ 

(13) 

where g u and gi are the degeneracies of the upper and 
lower states, A u i and X U £ are the Einstein A and wave- 
length for the transition, b is the Doppler parameter, 5x 
is the abundance relative to hydrogen for the element 
in question, and Ay t o — Ay/1 mag. The three num- 
bers given in square brackets are the numerical results 
for the lines [C n] 158 ^m, [O i] 63 ^m, and [O i] 145 
fim, respectively, using abundances 5c = 1-1 x 10~ 4 Z' 
and 6o = 5.0 x 10~ 4 Z' . This implies that, for the great 
majority of the cases I consider, optical depths effects 
will have at most a marginal effect on the cooling rate 
and can thus be neglected. For optically thin cooling 
at densities well below the critical density, the radiative 
cooling rate per H nucleus is simply 

Ai™ = Y] ^5 l E u , l e- E ^/ kBT C{k l ^nYi + h^ e n e ), 

(14) 

where the sum runs over the three upper states for the 
cooling lines (the 2 Py 2 state of C n and the 3 Pi and 3 Pq 
states of O i), g U: i and gg t i are the degeneracies of the 
upper states and the corresponding ground states, E u ^ 
is the energy of the upper state relative to ground, Si is 
the abundance of the relevant species, n e is the free elec- 
tron density, and fc.;_H and fci_ e are the rate coefficients 
for collisional de-excitation of the level by H and by free 
electrons, respectively. The clumping factor C appears 
for the same reason as for H2 formation on dust grains: 
these are collisional processes whose rates vary as the 
square of the local volume density. Obviously once H2 
becomes dominant over H 1 one should consider collisions 
with H2 as well, but since conversion to H2 only happens 
long after the gas has reached thermal equilibrium, I ig- 
nore this complication. 

I take the free electron density to be n e = 6cnn, 
which assumes that singly ionized carbon is the dominant 
source of free electrons; with this choice, excitation by H 
generally dominates. At low column densities, free elec- 
trons produced by ionization of hydrogen by soft x-rays 
outnumber those coming from carbon, but in regions of 

column density A^h ^ 1 20 cm" 2 this sourc e of electrons 
becomes subdominant (|Wolfire et al.l 120031) . Since most 
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Fig. 1. — Time evolution of the temperature (black solid lines) 
and H2 fraction (blue dashed lines) for a cloud of initial tempera- 
ture T = 1000 K, density no = 30 (clumping factor C\ = 1), and 
extinction Ay = 2 mag, initially composed of pure H I. Gray dot- 
ted vertical lines indicate times of t = tg and t = 10tg. Different 
lines are for metallicities of logZ' = —4 — in steps of 0.5, as 
indicated. 

of the cases with which I am concerned are in this regime, 
I adopt the carbon-dominated limit. 

Combining all heating and cooling processes, the tem- 
perature evolution obeys 

3 d 

^k B -T = r PEfi Z'e-^ N « + CCR9CR 
- V 9 —(h-K + Sch^E^e-^/^Cnv. (15) 

3. MODEL RESULTS 

3.1. Fiducial Parameters 

In order to survey parameter space, I consider a grid 
of model clouds of density Sh = 10° — 10 3 cm -3 in steps 
of 0.1 dex, extinction Ay = 10~ 2 — 10 1 mag in steps 
of 0.1 dex, and metallicity Z' = 10~ 4 — 10° in steps of 
0.05 dex, using the fiducial values for radiation and cos- 
mic rays specified above. All model clouds start with 
/h 2 = and T = 1000 K0 I integrate each model for 
15 Gyr and record the properties at t = tg, t = lOtg, 
and at t = 15 Gyr, with the latter representing the equi- 
librium state attained after long times; obviously this 
is longer than the age of the Universe, but these mod- 
els are useful to help build intuition for the importance 
of non-equilibirum effects. From the temperatures pro- 
duced in the models, I compute the Bonnor-Ebert mass 

Mbe = 1.18c 3 /-\/G 3 ^hTOhw, where c s — ^JksT / [im-n is 
the sound speed and ^ w 1.3 is the mean mass per par- 
ticle (as opposed to the mean mass per H nucleus /Up- 
values of Mbe should serve a rough proxies for where 
star formation can occur, with values much larger than 
the mass of any star indicating little star formation, and 
small values indicating star formation. 

Figure [TJ shows the thermal and chemical evolution of 
some example clouds drawn from the model grid. The 
figure is consistent with the qualitative timescale esti- 
mates above: gas at an initially high, non-equilibrium 

2 Obviously some of these initial conditions (e.g. clouds with 
Wh = 10 3 cm -3 and T = 1000 K) are unlikely to be found in real 
galaxies, but such models constitute a small fraction of the param- 
eter space, and the goal of this study is to perform a broad survey 
rather than trying to focus on particular assumed sets of "reason- 
able" initial conditions in galaxies where the actual conditions are 
poorly determined. 
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Fig. 2. — H2 fraction versus metallicity for model clouds, 
with blue indicating star-forming clouds, red indicating non-star- 
forming clouds, and green indicating intermediate values. For each 
model in a grid covering the parameter range log Z' = —4 — 0, 
log no = — 3, log [Ay /mag] = —2 — 1, I compute the Bonnor- 
Ebert mass Mbe and the H2 fraction at times t = iff (first panel), 
t = lOtff (second panel), and t = 15 Gyr (third panel). Blue pix- 
els indicate the locus in logZ', /h 2 corresponding to models for 
which Mbe < 100 Mq, indicating the star formation in the cloud 
is likely. Red pixels indicate the locus of models with Mbe > 1000 
Mq indicating star formation is unlikely; green pixels show inter- 
mediate values of A/be- F° r each color, brightness indicates what 
fraction of the models at that value of Z' have a given H2 fraction, 
with white representing none of the models, and solid blue, red, or 
green representing all of them. 

temperature will cool to a thermal equilibrium temper- 
ature of order 10 K and thus proceed to star formation 
in less than a free-fall time, even for metallicities as low 
as log Z' ~ —4. On the other hand, at metallicities be- 
low log Z' = — 1 the gas will be less than half converted 
to molecules at one free-fall time, and at metallicities of 
logZ' = —2 or less the gas will not reach 50% molecu- 
lar until more than lOtg. This result is consistent with 
numerical experiments in full cosmological simulations 
which show that equilibrium models of the H2 fraction 
begin to fail due to non-equilibirum effects at m etallicitcs 
below logZ' ss -2 (jKrumholz fc Gnediril20"Il . 

Figure[2]shows the H2 fraction as a function of metallic- 
ity for star-forming, non-star-forming, and intermediate 
clouds at t = tg, t = lOiff, and t = 15 Gyr (long enough 
that nearly all models have reached chemical equilib- 
rium). In clouds that are very old and thus have reached 
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Fig. 3. — Same as Figure [2] but computed for a UV radiation 
field 30 times stronger than the fiducial value. 

equilibrium, the figure shows that star-forming clouds 
(those with low -Mbe> indicated in blue) lie almost ex- 
clusively at high H2 fractions, and non-star-forming ones 
(those with high Mbe, indicated in red) almost exclu- 
sively at low H2 fractions, consistent with earlier work 
indicating a strong correlation b etween equilibrium gas 
temperature and chem ical state (jKrumholz et al.l 1201 It 
IGlover fc Clarkll2012al ). Out of equilibrium, the figure 
indicates that the correlation between low Mbe and high 
molecular fraction continues to hold for high metallic- 
ities. At lower metallicities, however, all models are 
displaced to smaller H2 fractions, and at metallicities 

log Z' ^ — 2 star-forming clouds are likely to have H2 
fractions well below unity even at t = 10tg. The impli- 
cation is that star-formation will be complete before the 
gas is significantly converted into H2 . The precise transi- 
tion metallicity below which equilibrium is not achieved 
will depend on the value of the clumping factor and the 
timescale for which star-forming clouds typically survive. 

3.2. Sensitivity to Parameter Choices 

To determine the sensitivity of these results to the 
choices of radiation field and cosmic ray intensity, I also 
compute the grid with a radiation field increased by a 
factor of 30 compared to the Milky Way value (Fpe.o = 
1.2 x 10~ 25 erg s" 1 and Cdi ss ,o = 1.5 x 10~ 10 s" 1 ) and 
with a cosmic ray flux increased by a factor of 10 to match 
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t = %, Cosmic rays x 10 



t = 10%, Cosmic rays x 10 



logZ' 



t= 10fff,z = 5 



log X 



log X 




Fig. 4. — Same as Figure [2] but computed for cosmic ray ion- 
ization rate 10 times larger than the fiducial value. 

the observed diffuse cloud value (Ccr = 2 x 10 _17 Z' s _1 ). 
The results are shown in Figures [3] and 01 Comparison 
with Figure[2]clearly indicates that the qualitative results 
are not substantially altered. 

Finally, note that in the thermal evolution calculation 
I have neglected heating due to cosmic microwave back- 
ground (CMB) photons. These will impose a tempera- 
ture floor T = 2.73(1 + z) K, where z is the redshift. A 
priori one would not expect the CMB to become signifi- 
cant until very high redshifts. The temperature reached 
by C II cooling does not fall below ~ 20 K over most 
of the model grid, and the CMB temperature does not 
exceed this value until z > 6.3. To confirm this intu- 
ition, in Figures [S] and [B] I show the results of imposing a 
minimum temperature T = 2.73(1 + z) K on the temper- 
ature used to evaluate A/be, for z = 5 and z = 10. The 
changes in the results from Figure [5] are essentially invis- 
ible at z = 5. At z = 10, the higher CMB temperature 
raises the temperature in some models such that there 
are fewer models with small values of Mbe, and these 
cluster at even higher molecular fractions. Qualitatively, 
however, the results are the same as at lower z. 

4. DISCUSSION 

4.1. Observational Implications and Tests 

The disconnect in timescales between H2 formation 
and cooling has two major observable consequences, 




log X 

Fig. 5. — Same as Figurc(2] but computed using a minimum gas 
temperature equal to the CMB temperature T^mb = 2.73(1 + z) 
K at z = 5. 

which can be used as a test of the above calculations. 
The first of these is a drop in H2 fractions below the levels 
predicted by equilibrium models in star-forming clouds. 
Figure [7] shows the ratio of the H 2 fraction at t = tg and 
t = lOiff to the equilibrium H2 fraction for star- forming 
clouds in the model grid. Clearly we expect equilibrium 
models to provide good predictions for galaxies down to 
metallicity Z' ss 0.1 or even somewhat less. This is con- 
sistent with observations to date, which show that chem- 
ical equilibrium models provide ex cellent fits to observed 
H 2 to Hi ratios in the Milky Way (jKrumholz et al.ll2009t 
I Lee et all 120121 ) and even the Sma ll Magellanic Cloud 
(SMC; Z' » 0.2. IBolatto etalllMl . However, the Fig- 
ure indicates that at metallicities of \ogZ' = —3, the 
H2 fraction in a given star-forming cloud will be at most 
~ 10% of its equilibrium value, and could be less than 
1% of that value. 

Second, the onset of star formation before the gas has 
time to fully transform to H2 in low metallicity galaxies 
should manifest as a reduction in the H2 depletion time 
idcp-H 25 defined as the ratio of the H2 mass to the star 
formation rate. In Solar metallici ty, non-starburstin g lo- 
cal galaxies, £dc P -H 2 ~ 2 Gyr (Bi giel et al.l 120081) . al- 
though lower values are possible in starbursts. This value 
should be lower in low metallicity galaxies by a factor of 
the mean H2 fraction in cold, star-forming clouds, since 
these clouds will only partially convert to H2 before form- 
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t = /,,, z = 10 




log T 




log X 



t= 15Gyr,z = 5 



Fig. 6. — Same as Figurc[2] but computed using a minimum gas 
temperature equal to the CMB temperature T^mb = 2.73(1 + z) 
K at z = 10. 




Fig. 7. — Ratio of H2 fraction at time t = tg (red) and t = 10% 
(blue) to the equilibrium value. For each model in a grid covering 
the parameter range log Z' = —4 — 0, log no =0 — 3, log Ay = 
—2 — 1, I compute the Bonnor-Ebert mass A/be H2 fraction at 
times t = tg and t = lOtfj. For this plot, I retain only models 
with Mbe < 100 Mq at t = tg, indicating that these clouds are 
likely to form stars. As in Figure[2] pixel brightness indicates what 
fraction of the models at a given Z' fall into that particular bin of 
/h 2 //h 2 e qi w ith white indicating no models and solid red or blue 
indicating all the models. 



Fig. 8. — Same as Figure[7] except that the y axis now represents 
the H2 depletion time tdep— H 2 = /h 2 (2 Gyr). 

ing stars and being destroyed by feedback. Figure [8] il- 
lustrates this effec t for clouds that live 1 a nd 10 free-fall 
times. Note that I Glover fc Clarkl ()2012b[ ) qualitatively 
suggested the existence of this effect, and Figure [5] rep- 
resents a quantitative extension of this prediction. 

Observational tests of these predictions are compli- 
cated by the fact that H2 is extremely difficult to ob- 
serve at low metallicities, because CO, the traditional 
H 2 proxy, ceas es to track H 2 at metallicities below a few 
tenth of Solar (iKrumholz et al.ll20lHlBolatto et alj|2011t 
Lerov et all 1201 U iShettv et all 120111: iNarav anan et al l 
20121 : iFeldmann et aLll20l3r Thus direct observational 
tests will require the detection of H2 by other means, 
such as dust or C 11 emission that is not associated with 
observed H I. While observationally challenging, surveys 
of this sort have alrea dy been comp leted in the closest 
galaxies like the SMC (jBolatto et alll201lD . and with the 
observational power provided by the Atacama Large Mil- 
limeter Telescope (ALMA) should begin to be possible 
in even lower metallicity nearby galaxies. In particular, 
850 /im observations are an excellent probe of dust and 
thus all gas including H2, because at 850 /jm dust is gen- 
erally optically thin, the emission is not very sensitive 
to dust temperature, and ALMA can achieve both high 
spatial resolution and excellent sensitivity. Prime tar- 
gets for such a campaign include IZwl8, SBS 0335-052 
(bot h Z' k, 0.02, and probably even lower dust me tallic- 
ities. Ilzotov et allll999tlHerrera-Camus et al.ll2012T ) . and 

Leo T (Z' < O.Ol. lSimon fc Gehal[2007l) . The ALMA ob- 
servations will have to be coupled with high resolution, 
high sensitivity H 1 maps to measure the atomic con- 
tent. Fortunately, su b-kpc resolution H 1 maps of IZwl8 
(jvan Zee et al.lflflQSl) and SBS0335-052 (lEkta et al.ll2009fl 
are already available in the literature. 

4.2. Implications for Simulations and Semi-Analytic 

Models 

These results have important theoretical implica- 
tions as well. Many galaxy simulation models al- 
low star formation only in regions where the gas has 
converted to H2; some of these models include non- 
equilibrium chemistry for H 2 form a tion and destruction 
dPelupessv fc Papadopoulod 120091: iGnedin et al l 120091: 
iGnedin & Kravtsovll2010HChristensen et al.ll2012D. while 
others assume equilibrium (|Fu et all 120101 : lLagos et all 
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[201lHKuhlen et al.ll20H Krumholz fc Deke]|[20T^ . The 
non-equilibrium models on average yield less H2 and thus 
less star formation at low metallicity, because often gas 
clouds are not able to build up significant H2 fractions 
before being dest royed by galactic shear or s imilar kine- 
matic processes ({Krumholz fc Gnedinl l201ll ) . However, 
if the relevant timescale is the cooling time and not the 
H2 formation time, and this effect should be far less sig- 
nificant. As a result, star formation should in fact oc- 
cur even in gas with low H2 fractions, provided that the 
equilibrium H2 fraction is high - it is the equilibrium 
H2 fraction and not the instantaneous one that corre- 
lates with gas temperature and thus is a good predictor 
of where star formation will occur. This suggests that, 
ironically, models in which the H2 is assumed to be in 
equilibrium, while they are less accurate in predicting 
the actual H2 fraction, may in fact be more accurate that 
the non-equilibrium models in predicting where star for- 
mation should occur. More generally, the calculations 
presented here suggest that star formation thresholds in 
simulations should be based on the instantaneous den- 
sity and extinction, which determine the temperature, 
and not on non-equilibrium chemical abundances. 

5. SUMMARY 

I explore under what conditions and for what physical 
reasons the observed correlation between star formation 
and molecular gas in the ISM is likely to break down. 
I show that the breakdown occurs at metallicities be- 
low a few percent of Solar, and that the physical mech- 
anism for this breakdown is a disconnect between the 
thermal and chemical equilibration timescales. Carbon 
in the ISM is able to cool gas on a timescale shorter by 
a factor of several thousand than that required for dust 
grains to convert the H 1 to H2. As long as both the 
thermal and chemical equilibration timescales are short 
compared to cloud free-fall times, which is the case at 
Solar metallicity, this does not have any practical effect 



and non-equilibium chemistry is unimportant. However, 
both the thermal and chemical timescales scale linearly 
with the metallicity, while the free-fall time time does 
not. At metallicities below a few percent of Solar, the 
free-fall time becomes intermediate between the thermal 
and chemical timescales, and clouds cool and proceed to 
star formation before molecules form, breaking the H2- 
star formation connection. 

This result has three major implications, two obser- 
vational and one theoretical. The observational impli- 
cations are that the equilibrium chemistry models that 
perform extremely well in the Milky Way and the SMC 
should begin to overpredict H2 abundances in very low 
metallicity galaxies, and that star formation should oc- 
cur in atomic-dominated regions of such galaxies as well, 
leading to a lower H2 depletion time. These predictions 
are not trivial to check, given the difficulty of measur- 
ing H2 in low metallicity environments, but combining 
high resolution dust and H I maps to infer the presence 
of H2 constitutes a viable strategy. The theoretical im- 
plication is that galaxy evolution simulations and semi- 
analytic models that link star formation to the chemi- 
cal state of the gas, and that treat that chemistry using 
non-equilibrium models, are likely to underpredict star 
formation rates in circumstances where the gas should 
reach thermal but not chemical equilibrium. It is the 
former that matters for star formation, not the latter. 
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